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ABSTRACT 



Aims. The first star forming galaxies in the early universe should be copious Lyor emitters, and may play a significant role in ionizing 
the intergalactic medium (IGM). It has been proposed that the luminosity function of Lya emitting galaxies beyond z~6 may be used 
to constrain the neutral fraction of the IGM during this epoch. In this work we report on a search for Lyo- emitters at redshift 8.8. 
Methods. We performed a narrow band imaging programme using ISAAC at the ESO VLT. Seven fields, covering a total area of 31 
sq. arcmin and for which optical and broad band infra-red images have been obtained in the GOODS survey, were imaged to a limiting 
flux (respectively luminosity) of ~ 1.3 x 10"'^ ergs s"' cm"- (respectively ~ 1.3 x 10*^^ ergs s"') in a narrow band filter centered in a 
region of low OH sky emission at 1.19yum. Candidate Lymana emitters are objects that are detected in the ISAAC NB images and 
undetected in the visible broad band images. 

Results. No z = 8.8 Lyo- emitting galaxies were detected to a limit approaching recent estimates of the luminosity function at z ~ 6. 
Our results do suggest, however, that detections or substantial constraints could be achieved by this method in the near future with 
larger field instruments planned for various telescopes. 

Key words. Galaxies: luminosity function - (Cosmology): early Univserse 



1. Introduction 

The reionisation of the universe is intricately linked to the for- 
mation of the first objects. While observations of high redshift 
QSOs indicate that the end of reionisation occurred around a red- 
shift of 6 (Fan et al. 2002), the measurement of the optical depth 
to electron scattering in the Wilkinson Microwave Anisotropy 
Probe (WMAP) three year results indicates that the u niverse was 
partially ionised up to s ignificantly higher redshifts dPage et alj 
12006 ^. 'Sperge l et al]l200 6). Hence, the current data point to an 
extended period of reionisation, ending at a redshift of 6 and ex- 
tending to redshifts well above 10. 

The objects that cause reionisation need not be the same 
throughout this period. A plausible model that matches the 
WMAP three year data is a model in which Population II 
and Population III o bjects contribute equally to reionisation 
dWvithe & Cenll2006h . In this model, the ionised fraction stays 
approximately constant over the redshift interval z - 7 - 12. 
However, other models, leading to quite different ionised frac- 
tions in this redshift interval, are also possible, so observational 
constraints on the fraction of ionised gas are now quite impor- 
tant. 

The first star forming objects in the universe should be co- 
pious emitters in Lyor and may play a significant role in ionis- 
ing the intergalactic medium (IGM). It has been proposed that 
the luminosity function of Lya emitters (LAEs) may be used to 
constrain the neutral fr action of the IGM when th e neutral frac- 
tion is as high as 30% (iMalhotra & Rhoadsll2004 . since the red 
wing of the Gunn-Peterson trough should attenuate the Lya line. 
As a comparison, the Gunn-Peterson trough measured blueward 



of the Lyg line saturates when the neutral fraction reaches 1 % 
(iFan et alJl20"02l) . 

Current results on the luminosity function of LAEat z = 
6.5 are contrad ictory. While Malhotra & Rhoadi (|2004 and 
iTaniguchi et al.l (l2005l) find no evidence for evolution in the LAE 
luminosity function between z - 5.7 and and z = 6.5, thus sug- 
gesting that the univers e was still largely reionised at z = 6.5, 
iKash ikawa et al. (2006) find significant evolution. 

Narrow-band, wide-field imaging surveys are an efficient 
and effective method of discovering LAEs. The narrow band 
filters are designed to cover regions where the night sky 
background is low. This naturally leads to surveys occur- 
ring at discrete redshifts, such as z = 3.4, corresponding to 
5390 A (iCowie & Hulll998l): z ^ 4.5, corresponding to 6740 A 



* Based on observations collected at the European Southern 
Observatory, Chile, Programs 072.A-0454, 076.A-0775. 



(iHu et al.lll998l: IMalhotra & Rho ads 2002); z^5.7, correspond- 
ing to 8150 A (Rhoads et al. 2003; Aiiki et al. 2006|; iHu et al.l 
2004; Shimasaku et al. 2006; Westra et al 2006) and z = 6.5 
corresponding to 9200 aIHu et al.l20023gi;lKodaira et al.ll2003l; 
ICubv et aI.ll200llTaniguchi et al.ll2005h . 

To detect Lyman alpha emitters (LAE) at still higher red- 
shifts, one has to move into the near IR. The next two windows 
that are relatively free of bright OH lines occur at L06//m and 
1.19;um, and correspond to LA Es at z = 7.7 and z = 8 .8, re- 
spectively. In a pioneering study, IWillis & CourbinI (|2005|) used 
ISAAC on Antu (VLT-UTl) to image a single field in the Hubble 
Deep Field South with a narrow band filter centred 1.19 fim. In 
a volume of 340 h^^ Mpc^ they find no LAEs brighter than the 

survey limit of IC*^ /i"^ ergs s"'. 

In this paper, we complement the work of IWilhs & CourbinI 
(l2005l) bv using the same instrument and filter to image a wider 
region of the sky to a brighter survey limit. In sections 2 and 
3, we describe the new naiTow band IR imaging data obtained 



2 



J.-G. Cuby et al.: A narrow-band search for Isya emitting galaxies at z = 8.8 



in this survey and how we combine these data with broad band 
optical and IR data from the GOODS survey to search for LAEs 
at z = 8.8. In section 4, we describe spectroscopic follow-up of 
our most promising candidate, and in section 5 we estimate the 
likelihood of detecting supernovae and solar system objects as 
false candidates. In section 6, we place limits on the z = 8.8 
luminosity function, and we discuss the possibility of extending 
this study over a wider area. Throughout this paper we use a flat 
ACDM model with Q.m = 0.3 and Hq = 70kms ' to compute 
distances and volumes. 



2. ISAAC observations and data reduction 

2.1. ISAAC narrow band observations 

We used the short wavelength (SW) arm of ISAAC 
dMoorwood et al.llI998l) on Antu (VLT-UTl) to image selected 
fields with the ISAAC narrow band filter centred at 1.19 micron. 
The detector in the S W arm of ISAAC is a Hawaii HgCdTe array, 
the pixel scale is 0'.'148 and the field-of-view is 2.'5 by 2'5. 

All the data were taken in service mode over a period last- 
ing several months. The observations started on 2004 November 
1 3 and finished on 2005 March 2 1 and occurred on 20 different 
nights. 

To ensure background limited performance, each exposure 
lasted approximately 20 minutes. After each exposure, the tele- 
scope was off'set by approximately 20" and a new exposure 
taken. This sequence was repeated while conditions were good 
enough and typically 3 to 6 exposures were taken on any given 
night. 

We targeted seven ISAAC fields within the Chandra Deep 
Field South (CDF-S) field covered by the Great Observatories 
Origins Deep Survey (GOODS). A list of the fields and the cor- 
responding exposure times and limiting magnitudes are given in 
Table[T] The total area of the survey (defined as being the sum of 
the areas that were continuously imaged) was 3 1 square arc min- 
utes. This area has been extensively studied at multiple wave- 
lengths - from the X-ray to the mid-IR. For the purposes of this 
study we make use of publicly available near-IR (VLT/ISAAC) 
(version 1.5) and optical (HST/ACS) data (version 1.0). 



Table 1. The fields observed and their coordinates, the expo- 
sure times, the magnitude at which the narrow band data are 
50% complete and the image quality of the stacked images. All 
magnitudes are in Vega and have been coiTected to 4" diame- 
ter apertures. The 90% completeness limit is approximately 0.5 
magnitude brighter. The field names correspond to those used at 
|http ://w w w. eso ■ org/science/goods/| 



Field 


RA 


DEC 


Exp. 


Lim. 


IQ 




(2000) 


(2000) 


time 


mag. 


FWHM 








(s) 


(Vega) 


(") 


F09 


03:32:33.2 


-27:44:55 


10800 


23.4 


0.55 


FU 


03:32:10.6 


-27:44:55 


18000 


23.2 


0.46 


F14 


03:32:38.8 


-27:47:25 


21600 


23.2 


0.55 


F15 


03:32:27.5 


-27:47:25 


6000 


22.4 


0.69 


F16 


03:32:16.2 


-27:47:25 


26400 


23.2 


0.53 


F20 


03:32:33.2 


-27:49:55 


18000 


23.3 


0.59 


F21 


03:32:21.9 


-27:49:55 


20400 


23.3 


0.58 



2.2. ISAAC narrow band reduction 

Apart from the removal of the background, the methods used 
to reduce the ISAAC narrow band data were standard. A dark 
frame was first subtracted and a twilight flatfield applied. The 
removal of the background was non-standard and had to be done 
in two steps. Usually, when processing data in the near IR, the 
background of any particular frame is estimated by averaging the 
frames that are taken immediately before and after the frame in 
question. To maximize the signal-to-noise ratio, one should use 
as many frames as possible, but this is limited either by the num- 
ber of frames that are available or by the timescale over which 
the background varies. In these data, the large scale component 
of the background was varying with a timescale that was similar 
to the length of the exposure. 

The reason for the variation is that, within the field of view of 
ISAAC, the central wavelength of the narrow-band filter trans- 
mission curve varies. In the upper part of the field, light from 
OH lines pass through the shifted filter band pass. Since the in- 
tensity of OH lines can vary considerably over a period of a few 
minutes and since exposures were relatively long (20 minutes), 
there was considerable variation in the large scale background, 
especially in the upper part of the array. 

In addition to the relative rapid changes in the background, 
sometimes, there were as few as three frames taken consecu- 
tively, so, in principle, one had only two frames to use as sky. 

In order to remove the background adequately and to maxi- 
mize the signal-to-noise ratio in the background subtracted im- 
ages, we first removed the large scale shape of the background 
in individual frames by applying a median filter with a box size 
of 200 pixels to each image and by subtracting the result. Once 
the large scale shape is removed, we then collect all the data 
for a given field, even if it was taken on different nights, and 
follow the standard procedure of background subtraction that is 
descibed above. This unusual approach only works if the small 
scale variations (including pixel-to-pixel variations) are stable 
over the time the data are taken. This appears to be the case for 
these data as the pixel-to-pixel noise in the reduced data is within 
a few percent of the theoretical expectation. 

We used the XDIMSUM package in IRAlQ to do the back- 
ground subtraction as it allows object and bad pixel mask- 
ing when computing the background frames. Typically, a back- 
ground frame is made up of 15 to 20 frames. 

The alternative of subtracting the backgound first and then 
removing any large scale residuals with a median filter proved 
to be unsatisfactory, as objects, variable pixels and cosmic rays 
cannot be properly identified by the XDIMSUM software when 
the background variation between frames is as large as it was in 
these data. 

The background subtracted images were then corrected for 
field distortion, using the field distortion coefficients available 
from the ISAAC web pagefl registered, combined and astro- 
metrically aligned to the broad band ISAAC and ACS images. 
Because some of these steps involve interpolation, the pixel-to- 
pixel noise in the fully reduced images is a few percent less than 
the value one would derive from the mean level of the back- 
ground. 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement with 
the National Science Foundation. 

^ http://www.eso.org/instruments/isaac/ 
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3. Analysis 

3.1. Catalogue generation and detection limits 

We used the SExtractor (Bertin & Arnouts,1996) software in sin- 
gle image mode to detect and measure objects in the ISAAC 
narrow band images and in the ISAAC and ACS broad band im- 
ages. We created one catalogue per filter per field and we then 
used our own software to match the catalogues in sky coordi- 
nates. The matching radius was set to between 0.2 and 0.4 arc 
seconds, depending on the field. 

The magnitudes were first computed for apertures of 8 pixels 
in diameter and were then corrected to apertures of 4" in diam- 
eter by applying an aperture correction that was computed from 
bright unsaturated stars. 

The limiting magnitude of the naiTow band data can be esti- 
mated by comparing the number of detections in the narrow band 
image with the number of detections in the much deeper ISAAC 
J-band image. Since the ISAAC image is much deeper and since 
J - NBl 190 ~ for most objects, objects detected in the narrow 
band filter should also be detected in J. By plotting the ratio of 
the number of objects seen in both the narrow band image and 
the J-band image over the number of objects seen in the J-band 
image, one can get an estimate of the completeness. An example 
for one of the fields is shown in Fig.[T] The completeness limit is 
defined as the magnitude at which the narrow band data is 50% 
complete and is estimated by fitting the function 



f(m) = (exp(^^^— p-l-H 1 



(1) 



to the completeness fraction as a function of magnitude. The 
parameters and are the 50% completeness limit and a 
measure of how quickly incompleteness sets in, respectively. 
The 50% completeness limits of the survey fields are listed in 
Table [T] For an object with magnitude nic, the signal-to-noise 
ratio over an 8 pixel diameter aperture is approximately 6. In 
computing this signal-to-noise ratio, the coiTelation in the noise 
between pixels has been taken into account. By creating an artif- 
ical image that has the same pixel-to-pixel noise as the data, one 
finds no false detections to the same detection threshold. 

Objects that are detected only in the narrow band image are 
either line emitting objects with large equivalent widths, ex- 
tremely red objects or false detections. In Fig. [1] these objects 
are plotted in dot-dashed line. Their number does not become 
significant until one passes the 50% completeness limit. Beyond 
this limit, most objects that are only seen in the narrow band 
image are probably false. So we choose this limit as the limit 
of our survey. With the exception of field F15, which is signifi- 
cantly shallower than the other fields, the average limiting mag- 
nitude is 23.3, corresponding to a flux (resp. luminosity) limit of 
~ 1.3 X 10^'^ergss"' cm"^ (resp. ~ 1.3 x lO'^-' ergs s"'). The flux 
and luminosity limits do not account for partial line transmission 
effects. 

3.2. Candidate selection 

Lyman alpha emitting candidates are selected using the follow- 
ing criteria. They should be 

- visible in narrow band ISAAC image and 

- not visible in the optical (F435W, F606W, F775W, and 
F850LP) ACS images 

We note that these criteria do not make use of the IR data: 
no criterion was set a priori for the equivalent width of Lya at 




Magnitude 



Fig. 1. Panel (a): A comparison of the number of objects de- 
tected in the J-band image (solid line), in both the narrow band 
and J band images (dashed line), and in the narrow band image 
alone (dot-dashed line). Panel (b): A histogram of the complete- 
ness function of objects detected in the narrow band filter. The 
histogram corresponds to the ratio between the dashed and solid 
lines in the upper panel. The dotted line is a fit to the histogram 
using equation[T] 



high redshift, and the IR data were used a posteriori to remove 
low redshift objects (see below). This step led to an initial list 
of ten candidates which we then visually inspected in all bands, 
including the ISAAC J, H and K band data. Down to the limiting 
magnitudes listed in Table [1] all candidates were found to be 
either 

- ghosts related to bright stars (two objects) 

- in areas of low signal-to-noise that are near the edges of the 
combined narrow band images and are, therefore, outside our 
survey area (one object), or 

- very red objects not detected in the optical bands with sig- 
nificant flux in the J, H and K ISAAC bands (five objects) 

In the later category can be found Extremely Red Objects 
(EROs) (R - K > 5) or T-dwarfs which exhibit even redder col- 
ors (I-J > 6). These objects can therefore pass our initial selec- 
tion criterion if they are detected in the NB image. ER Os have 
increasingly red colors from NB to K (J - K > 1.8) tCimatti et al.l 
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l2003h . whereas T-dwarfs have approximately flat J, H and K col- 
ors and NB (1.19jt/m) - J ~ 1 due to a strong water absorption 
band at 1.19 //m. We then inspected all initial candidates and 
further eliminated those which had broad band IR colors consis- 
tent with one or the other of these very red objects. In practice, 
only objects consistent with ERO colors were found. Note that 
the data points of the dot-dashed line of figure 1 correspond to 
the candidates before this second step was performed. 

After this careful refinement of the candidate list, we were 
left with no objects that were brighter than the 50% detection 
limit, and two candidates fainter than this limit, both only vis- 
ible in the NB1190 image. In view of the limiting magnitudes 
of the J and NB images, this corresponds to 7 - A^B1190 ^ 1.5, 
which in turn corresponds to objects with observed (resp. rest- 
frame) equivalent widths EWohs. ^ 300 A (resp. EW^si ~ 30 A) 
which are acceptably large values for a very high redshift star- 
burst galaxy. 

Even though the magnitude of the two candidate were fainter 
than the 50% completeness limit, we decided to observe the most 
promising one (NBl 190=23.9 and S/N=3 in the F21 field) with 
SINFONI with the aim of verifying that the object was real. 

4. SINFONI observations 

We used SINFONI on Yepun (VLT-UT4) on the night of 2005 
October 31 to observe the most promising candidate in the 
survey. SINFONI is a n ear infrared integral field spectrograph 
dEisenhauer et al.ll2003h that can be fed with an adaptive optics 
(AO) corrected beam. Since there were no bright stars near to 
the candidate, we used SINFONI without AO. We used the 8" 
field-of-view, resulting in a pixel scale of 0'.' 125 perpendicular to 
the image slicer. The image slicer splits the field-of-view into 32 
slices, so each slice is 0'.'25 seconds wide. The 32 slices were dis- 
persed with the J-band grating and imaged onto a Hawaii 2kx2k 
detector. The resulting spectra start at 10950 A and effectively 
end at the beginning of the strong telluric absorption feature that 
starts at ~ 13500 A. The resolution of the spectra is about 2000. 

Individual exposures lasted 900 seconds, which is long 
enough to be detector dark noise limited at 1.19 microns. A sin- 
gle observing sequence consists of four such exposures with the 
object dithered within the 8" field-of-view. At the end of the 
sequence, a nearby star, which was used to verify the pointing, 
was observed. This sequence was repeated 6 times and results in 
a total integration time of 21600s. 

The observing and reduction strategies were designed to 
maximise the signal-to-noise ratio of a point-like emission line 
object. The first step in the reduction process was to create a 2 
dimensional sky frame for each exposure. For any given expo- 
sure, the corresponding sky frame was made from the average of 
the 6 other frames. Thus, instead of reducing the optimal signal- 
to-noise ratio by a factor V2, which would be the case if the sky 
subtraction was done on the basis of image pairs, the signal-to- 
noise ratio is reduced by ~ 8%. The maximum number of frames 
that one can use for a sky is limited by instrument flexure. We 
found that 6 frames, taken over a period of two hours, was the 
best compromise between maximising the signal-to-noise ratio 
and minimising the residuals caused by flexure. 

The sky-subtracted data were then flat fielded and sectioned 
into 32 images, one for each slice. The 2 dimensional sky sub- 
traction is not perfect because, as a function of wavelength, the 
relative intensity of the IR night sky varies with time. Night 
sky residuals were removed by fitting and subtracting low order 
polynomials along the spatial direction of the slices. 



The resulting 2 dimensional sky-subtracted data were then 
registered and averaged to produce a 2 dimensional co-added 
frame. A one square arc second region was then extracted from 
the co-added frame to produce a 1 dimensional spectrum that 
was then wavelength and flux calibrated. 

After six hours of integration, no emission line was detected. 
Our 5 sigma detection limit over a one square arc-second aper- 
ture is 1.5x10"'^ erg s"' cm"^ pixel"', or 3.0x 10"'** erg s"' cm"^ 
if we integrate over 6 Angstroms (4 pixels). From the magnitude 
of the candidate, we expected a flux of 7.0 x 10""* erg s"' cm"^. 
We conclude that our best candidate was either spurious or a 
transient object (see next section). 

5. Contamination from transient objects 

Since the narrow band and broad band data were taken at dif- 
ferent times, it is possible for transient objects to appear in one 
filter and not another. Our candidate selection is based on find- 
ing sources that can only be seen in the ISAAC narrow band im- 
ages. Therefore, transient objects that are brighter than the nar- 
row band detection limit will be considered as candidates if they 
are not visible in the broad band images. We consider two tran- 
sients that might contaminate the survey: supernovae and bodies 
in the solar system. 

A few of the fields were observed at a time when the angular 
motion of distant solar system bodies as seen from the Earth 
can be very low. In individual exposures, which typically last 
20 minutes, such objects may appear point-like. Since the total 
exposure times were long and since the observations were done 
in service, the data were taken during different nights. For field 
21, the time over which the data were taken spans 4 days; for 
field 16 it took 6 weeks to collect the data. This, together with 
the location of the fields, which are 45 degrees from the ecliptic, 
means that it is very unlikely that we would mistakenly identify a 
distant body in the solar system as a high redshift galaxy. When 
the data are combined, objects in the solar system will either 
be rejected when clipping for cosmic rays or appear as trailed 
irregular-looking sources. In either case, it is sufficient to split 
the data into subsets that are based on the date of the exposure to 
test if a candidate is moving or stationary. We have done this for 
all our initial candidates and no moving targets were identified. 

At the flux limit of the survey, distant supernovae can 
be visible for several weeks, so supernovae are a potential 
source of contamination. We computed the expected number of 
Type la and Type II supernovae that would be visible in our 
nai Tow band ISAAC imag es by using the rates published in 
Ho pkins & BeacomI (|2006|) Even though Type II supernova are 
approximately 5 times more numerous than Type la supernova, 
they are, on average, 2.5 magnitudes fainter ( Filippenko 1997|). 
so Type la supernovae dominate the expected counts. Using a 
magnitude limit of 23.5 in the narrow band filter, we would need 
approximately 100 square arc minutes to detect one Type la su- 
pernova. This is only three times the area covered in our survey, 
so it is possible that there are Type la supernova that are visible 
to the limit of the survey. However, Type la supernova are usu- 
ally close to their host galaxies. Since most of these hosts would 
be detected in our broad band images, most Type la supernovae 
would not be selected as candidates. There are, however, a num- 
ber of supernova that are considered as hostless. I.e. to the limit 
at w hich imaging data is available, no host is visible. From Table 
4 in Lidman et al.l (12005^ one estimates that 1 in 10 high redshift 
Type la supernova can be considered as hostless. However, the 
GOODS ACS imaging data is considerably deeper than what is 
usually taken when following high redshift supernovae, so the 
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fraction of hostless Type la supernova is likely to be even less 
than 1 in 10. 

Although it is very unlikely that we would have selected a 
supernova of any type as a candidate in our survey, we can- 
not completely rule it out. In surveys that cover larger areas or 
use less extensive broad band imaging, contamination from tran- 
sients such as supernova or more unusual events, such as the one 
recently discovered in Bootes (Dawson et al. 2006), might occur, 
so strategies to avoid this contamination, such as taking the nar- 
row band data over several epochs separated by several months, 
should be considered. 

6. The z ~ 8 LAE luminosity function 

No Lya emitters at z ~ 8.8 have been found in 31 arcmin^ down 
to ~ 1.3 X 10~'^ ergs s~' cm ~^. Our results complement the work 
of lWilhs & CourbinI (l2005h who performed similar observations 
with the same filter and instrument, however with different sur- 
vey area and depth. 

From a large sample of 58 Lya emitters (LA Es) including 17 
spectroscopic confirmations, iKashikawa et al.l (,2006.) have es- 
tablished the luminosity function of z = 6.5 LAEs. Similarly, 
IShimasaku et al.l ([2006) have established the luminosity func- 
tion (LP) of z = 5.7 Lya emitters from a sample of 89 LAEs, 
of which 28 spectroscopica l ly confirmed. UnUke pre vious work 
dMalhotra & Rhoadsll2004l) . IKashikawa et al.1 (|2006|) show con- 
vincing evidence for evolution of the LP between z = 5.7 and z 
= 6.5. It is not clear whether this evolution can be attributed to 
the luminosity L*, the comoving number density O*, or the faint 
end slope a, or a combination of those, although it is suggested 
that L* is dimmed by a factor ~ 2 between these two redshifts. 

In figure|2]we represent the LAE luminosity functions at red- 
shifts 5.7, 6.5 and at z > 7 assuming arbitrarily another dimming 
in L* by a factor 2 between redshifts 6.5 and > 7. Also repre- 
sented on t his figure are the data po ints corresponding to this 
work and to ^ Willis & CourbinI (|2005|PI The two data points are 
consistent with the likely evolution of the luminosity function 
towards higher redshifts, but do not allow us to constrain this 
evolution due to the limited number of fields covered. We also 
indicate on figure |2] a point coiTesponding to an on-going pro- 
gramme that we are carrying out at CFHT with a wide field im- 
ager (WIRCAM) at redshift 7.7. When this programme is com- 
plete, we should either detect a few z = 7.7 galaxies if the lumi- 
nosity function has moderately evolved from z = 6.5, or, in case 
of no detections, infer a significant evolution of the LF. 

Other wide field imagers equipped with low OH NB fil- 
ters, either at 1.06 fim or 1.19 fim, are becoming avail- 
able on 4m or 8-lOm class telescopes. In 2007, the 
VLT will be e quipped with HAWK-I (see ICasaliet al] 
( |2006|) and |http://www.eso.org/ instruments/hawk i) , a ~ 8'x 
8' imager, and VISTA (see McPhers on eTaLl ('2006') and 
jhttp://www.vista. ac.uk), the new 4m IR telescope to be installed 

^ We note that all data used in figure |2| suffer from very large er- 
ror bars, not shown. One source of error comes from the way the 
comoving volume samp led by the NB imaging survey is computed. 
IWillis & CourbinI ( 120051) correctly compute the comoving volume as 
a function of the Lya luminosity (and line width), while other authors 
compute a fixed comoving volume from the full width at half maxi- 
mum of the NB filter. Such differences could significantly affect the 
luminosity function derived from number counts, and in particular the 
faint end slope of the LF. For consistency with other work, we have used 
the traditional way of computing the comovi ng volume when plotting 
our datapoint on figure |2l and we adapted the I Willis & Courbini ( |2005|) 
datapoint accordingly. 
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Fig. 2 . Luminosity fu n ction of LAEs at z - 5.7 (dashed 
line, Shimasaku et al. (2006*)) an d z — 6.5. (pla i n line , 
Kashikawa et al. (2006)). Square: IWilUs & CourbinI (l2005h . 
Circle: this work. Triangle: on-going survey at CFHT targeting z 
- 7.7 LAEs. Data points have not been corrected for complete- 
ness. Extrapolation of the LF to higher z (z > 7), assuming the 
same evolution from z=6.5 as between z=5.7 and z=6.5 is shown 
for illustration. 

at Paranal will cov er a field of view of ~ 45'x 45'. DAZLE 
(iHorton et al. I I20041) . a visitor instrument dedicated to searching 
high-z LAEs (7'x 7' field of view equipped with very narrow 
band filters at < 0. 1 % bandpass) will observe at the VLT in the 
coming months. 

These new instruments will be essential for investigating 
high-z LAEs by probing much larger volumes at similar or 
fainter detection limits. Different survey strategies can be con- 
templated with these instruments: either go very deep in a sin- 
gle field, or alternatively go brighter in several fields (pointings). 
Figure |3] shows the number of objects that can be detected us- 
ing two different survey strategies with HAWK-I at the VLT. 
Two luminosity functio ns are assumed: the z = 6.5 LF (L*) 
(IKashikawa et 31.1 12006*) and the ad' hoc z > 7 LF mentioned 
above (Lp. The two different strategies, for each of the two lu- 
minosity functions, compare the number of objects detected in 
a single 'deep' field and in 5 non-overlapping 'shallow' fields 
observed each for a 5 times shorter integration time. 

The plot indicates that the 'deep' strategy is naturally to be 
preferred for short integration times, while favoring area over 
depth becomes advantageous for very long integration times. 
Considering the very high uncertainties in the LF of z > 7 galax- 
ies, a single pointing appears preferable for NB imaging surveys 
up to 100 hrs of integration time or less for the instrumental char- 
acteristics considered in this example. 

7. Conclusion 

We have searched for z ~ 8.8 Lyor emitters with a NB filter cen- 
tered at 1 . 19 yum in ISAAC at the ESO VLT. We covered an area 
of 31 arcmin^ down to a line flux of ~ 1.3 x 10"'^ ergs s"' cm"^. 
No LAE was detected in this pilot survey, consistent with recent 
estimates of the luminosity function at z ~ 6. Similar observa- 
tions with wider field instruments are likely to rapidly unveil 
high-z LAEs at redshifts above 7, however in limited numbers 
and at high cost in observing time. Systematic surveys of large 
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Fig. 3. Number of detected LAEs versus time. The VLT/HAWK- 
I case is assumed. Two different strategies (one deep single 
pointing and 5 shallower pointings for the same total time on 
sky) are indicated, for t wo luminosity functions: the z - 6.6 LF 
dKashikawa et al.|[2006D (solid line, L*) and the same z > 7 LF 
as in figure 2 with dimmed by a factor two (dashed lines, Lp. 



numbers of very high z galaxies will have to await JWST and 
ELTs. 
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